Spin-Orbit Interactions in Bilayer Exciton- Condensate 

Ferromagnets 

Y.-P. Shimi'2 and A. H. MacDonald^ 

^Department of Physics, University of Texas at Austin, Austin, Texas 78712 
^Institute for Micro structural Sciences, 
National Research Council of Canada, Ottawa, Canada KIA ORG 

(Dated: May 4, 2009) 

Abstract 

Bilayer electron-hole systems with unequal electron and hole densities are expected to have ex- 
citon condensate ground states with spontaneous spin-polarization in both conduction and valence 
bands. In the absence of spin-orbit and electron-hole exchange interactions there is no coupling 
between the spin-orientations in the two quantum wells. In this article we show that Rashba spin- 
orbit interactions lead to unconventional magnetic anisotropics, whose strength we estimate, and 
to ordered states with unusual quasiparticle spectra. 
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I. INTRODUCTION 



An exciton is an elementary excitation of a semiconductor in which a conduction band 
electron and a valence band hole form a bound state. Like Cooper pairs in superconductors, 
excitons can condense under appropriate circumstances. The broken symmetry associated 
with exciton condensation is spontaneous phase coherence between conduction and valence 
band states. Excitonic condensation has been one of the most anticipated forms of fermion- 
pair condensation since first predicted in the early 1960's.->2i>^i^>^ 

Excitons are traditionally created by optically exciting valence band electrons to the 
conduction band. One key obstacle to the creation of an exciton condensate has been the 
finite life time of optically generated excitons, limited in most circumstances by optical 
recombination processes. It was proposed^"^ some time ago that the exciton life-time can 
be substantially increased by spatially separating electrons and holes in a bilayer configu- 
ration. Progress in layered semiconductor growth techniques has now made it feasible to 
realize bilayer electron-hole systems with great flexibility. Many experimental studies have 
attempted to detect condensation of optically generated excitons using photoluminescence 
measurements^'^i^iiiii^'i^ii^ in coupled quantum well structures. 

Systems in which electrons and holes are present in equilibrium can be achieved in suit- 
ably gated semiconductor bilayer systems^^'^^ and provide a simpler and more ideal realiza- 
tion of an electron-hole fluid. A degenerate bilayer electron-hole fluid undergoes excitonic 
condensation^ when the distance between the two layers is smaller than the average distance 
between particles in one layer. The interesting anomalous transport properties'^ associated 
with exciton condensate superfluidity can be studied only in equilibrium electron-hole sys- 
tems. The essential technical development necessary to perform transport measurements 
which probe exciton condensation is the perfection of techniques which enable separate 
electrical contacts to electron or hole layers.— Indeed magnetoexciton condensaton, which 
occurs in strong magnetic fields and involves electrons and holes that can be both in the 
conduction band, both in the valence band, or in separate bands, has already been realized 
and studied using transport.— i2Si2ii2^ Equilibrium bilayer electron-hole fluids can be created 
in weak magnetic fields by electrically generating an external potential difference between 
nearby quantum wells approximately equal to the energy gap of the host semiconductor. 
This paper is motivated by impressive recent progress in this direction.— One impor- 
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tant feature of electrically generated equilibrium electron-hole fluids is the ability to study 
the dependence of system properties on the difference between electron and hole densities. 
Because both conduction and valence bands are spin degenerate and condensation is favored 
by perfect Fermi surface nesting, it has been predicted that a finite population difference be- 
tween electrons and holes enriches the physics of excitonic condensates, leading in particular 
to spin-polarized phases.— 

In this paper, we consider a double-quantum-well structure with two GaAs quantum 
wells separated by an AlGaAs spacer. The larger band gap in AlGaAs compared to GaAs, 
leads to potential profiles in which AlGaAs acts as a barrier for both conduction band 
electrons and valence band holes. An equilibrium electron-hole fluid is created by applying 
a large electric field across the direction perpendicular to the layers so that the valence band 
maximum in one quantum well can move above than the conduction band minimum in the 
other. This procedure leads to spatially separated but strongly interacting electron and 
hole fluids. In addition to making electron-hole systems stable, the external field inevitably 
enhances structural inversion asymmetry and therefore introduces (Rashba) spin-orbit (SO) 
interactions^Ii^ in both quantum wells.— '^'^'^ Recent theoretical work— has shown that 
Rashba SO interactions lead to unconventional anisotropic electron-hole pairing in a bilayer 
system with equal densities of electrons and holes. The focus of this paper is on the role the 
Rashba SO interaction plays in the spin-polarized phase expected to accompany excitonic 
condensation whenever the electron and hole densities are unequal, which is referred to as 
the exciton-condensate ferromagnet state. 

We model the experimental system used to create the equilibrium electron-hole fluid 
by assuming that external charges outside of the bilayer system can be controlled to fix the 
electric fields E^,, E/j, on the left and right hand sides of the bilayer system. Treating E^ and 
E/j, rather than the electron and hole densities, as experimentally controllable parameters 
represents the real experimental condition more accurately. The system we consider is 
schematically shown in Fig. [H Two external layers {L and R) with charge density pi and 
pn controls the electric field E^, and E/j which determine the electron and hole densities in 
layer A and B. Electrons will be in layer A and holes will be layer B. We will use conduction- 
valence band picture rather than electron-hole picture to avoid confusion in calculating the 
total energy of the system. For layer A {B), we consider only the conduction (valence) 
band and neglect the valence (conduction) band. The charge density of layer A is given by 
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FIG. 1: (Color online) Cartoon depicting the bilayer system including the external charge distribu- 
tion which gives rise to the external electric field. Layer A is the electron layer and layer B is the 
hole layer. Layers L and R contain the external charge distribution. The overall charge is neutral 
so that pl + PA + Pb+Pr = ^- Pl and pn are purely external charge and lead to gate external 
fields Ei and E^. pA = —en^ is contributed by the conduction band electrons and ps = enh is 
the charge density contribution from valence band holes. The electric field outside of the system 
is zero due to the overall charge neutrality. 



Pa = (— e)ne. Layer B also contains background positive charge with density hq so that 
fully occupied valence band is charge neutral. The charge density of layer B is then Pb = 
eiriQ—n^) = erih- The overall charge neutrality is always preserved, i.e., Pl+Pa+Pb+Pr = 0. 
Once we fix the external charge densities pl and pr, the density difference between electrons 
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and holes An = Ue — Uh is determined by the overall charge neutrality requirement, but He 
and rih are individually determined only by requiring a constant chemical potential through 
the entire system. When and E/j are identical, Pl + Pb. = and therefore = Uh, while 
if El -Er ^ 0, then + 7^ and An = -{pA + Pi?)/e = {pi + P/?)/e. Therefore, 
by controlling the external charge densities pi and pr or equivalently, the external electric 
fields El and E/j, it is possible to induce density-polarized {rie 7^ n^) electron-hole bilayer 
systems. The ability to tune the difference between electron and hole populations is therefore 
a natural feature of electrically generated equilibrium electron-hole systems. This population 
polarization leads to interesting new physics that is our main interest. 

This paper is organized as follows. In Sec. [ITl we briefly review the mean-field theory of an 
exciton condensate for the case of spinless electrons and holes. In Sec. IIIH we restore the spin 
degree of freedom and study the spin-polarized states which occur when the electron and 
hole densities are unequal. We find that as the electron-hole population difference increases, 
the ground state progresses from one with pairing in two spin-channels, to a state with 
pairing in only one spin channel, and finally to a state with no pairing. States with finite 
population polarization generally have spin-dependent pairing. When spin-orbit interactions 
and electron-hole exchange interactions are neglected the ground state is invariant under ar- 
bitrary spin-rotations in conduction and valence bands. Each quantum well is spin-polarized 
but there is no coupling which favors any particular relative spin orientation. In Sec. IIVI we 
turn to a discussion of the role played by the Rashba spin-orbit interaction in the exciton 
condensate ferromagnet. We find that spin-orbit interactions can favor spin-alignment in 
a particular direction depending on the electron and hole densities and the strength of the 
SO interaction. In Sec. |V]we briefly summarize our findings and discuss some experimental 
implications of our work. 



II. SPINLESS EXCITON CONDENSATION 



In this section we establish our notation and introduce the basic equations which appear 
in the mean-field theory of electron-hole condensates by studying a system with spinless 
conduction and valence bands localized in spatially separated quantum wells. One of our 
main objectives here is to explain how we deal with the potentially confusing way in which 
simple electrostatic effects are combined with pairing in systems with electron and hole 
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densities which are spatially separated and controlled externally. As it turns out, the spin 
degree of freedom plays a minor role for unpolarized {ue = rth) systems, so the results 
in this section are applicable with minor revisions to electron-hole bilayers with balanced 
populations. 

We treat the charges in layer L and layer R as well as the background positive charge in 
layer B as external charges. The total Hamiltonian of a spinless bilayer system that has one 
conduction band and one valence band is 

H = E^il + (42 + 4c,k + ^ E 4 k Ca'k'+q C,k-q , (1) 

ak kk'q 

aa' 

where the band index a = c for the conduction band c in layer A, and a = v for the valence 
band in layer B. In Eq. ([T]) E'^^ is the electrostatic energy that comes from the Coulomb 
interaction between the external charges, and V^^^ is the electrostatic potential in band a 
due to external charges, k is the two-dimensional wave vector, f2 is area of each layer, and 
yaa' ^^iq Fourier transform of the Coulomb interaction between electrons in band a and 
band a' given by 

27re2 

V^\c^) = y^^(q) = ^ , (2) 
\/-(q) = \/-(q) = -^e-"'' , (3) 

where e is the dielectric constant. In bilayer systems, unlike single-layer systems, the Hartree 
electrostatic energy must be included explicitly in the mean-field theory since individual 
layers are not necessarily charge neutral. Only the overall charge neutrality is preserved. 
The bare band energies e^^ 's ^^^^ assumed to be parabolic for both conduction and valence 
bands, 

^(0) - ^ ^(0) - _ E (4) 

'^"^ ~ 2m, ' ^^i^ " 2m, ' ' 

with the fundamental band gap Eg. V^^^ does not depend on k and only depends on 

the 2;-coordinate of the layers. We neglect the small overlap between electron and hole 

wavefunctions so that the interaction V conserves band indices and that therefore there is 

no electron-hole exchange interactions. Allowing for spontaneous inter-band coherence, the 

mean-field Hartree-Fock Hamiltonian is given by 

Hmf = Eo + J2 (42 + <k + Vf 4cak - (Ak4c.k + h.c) , (5) 

ak k 



where 



ak 



(6) 



k' 

Ak = ^5^F-(k-k')p™(k') , 

k' 

P^^'i^) = (cLkCak)* = (^|cl,kCak|vl>) . 



(7) 



(8) 



(9) 



Note that the Hamiltonian itself depends on the state \1/ and must be solved for self- 
consistently. Here e^^ is the intra-band exchange field, Ak is the inter-band exchange field 
induced by inter-band coherence, and p""'(k) is the density matrix. These potentials, which 
must be also determined self-consistently, are discussed at greater length below. Ees is the 
electrostatic energy and V^^^ is the total electrostatic potential at layer a, due to the whole 
charge distributions in the system which includes the Hartree potential from the electrons 
and holes. The last two terms in Eq account for the double counting of the intra- and 
inter-layer exchange interactions in the mean-field Hamiltonian and must be included when 
we calculate the total energy. 

The intra-band exchange [Eq. ([7])] potential for the valence band appears to diverge 
because p™(k') 1 deep in the valence band. This problem is solved by noticing that 
the experimentally measured valence band effective mass m* already includes all interaction 
effects present in the state with a full valence band and an empty conduction band, a 
reference state we will define as |$o)- I'^o) is the solution of the problem with Pl=Pr=0. 



Thus 



hp 



4k^ + Ck[<^>o] 



(10) 



2m* 



he 



1 J]l--(k-k')pr(k') 



2m, 




(11) 
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where po is the density matrix when the system is in state |$o), 



, 1 a = V, a' = V , , 

pT{k)={ . (12) 

otherwise 



It follows that 



42 + ell = - - ^ E ^^^(k - k') (P"'(k') - pr(k')) , (13) 



which remains finite. The total energy of the state |\E') is measured with respect to the 
energy E^^l of the state |$o)- 

Our calculation requires a self-consistent evaluation of the difference between the density- 
matrix and the density matrix in the absence of carriers: 

p»'^'(k)=p-'(k)-pf (k) . (14) 

With this definition Eq. ([5]) becomes 

Hmf - El^l = Eo + Y, (4cak - pr(k)) - Yl (^k4c.k + h.c.) , (15) 



ak 



where 



^0 = E,s-EPs-Y^a^'fr^(^) 

ak 

-lY^ain^) + (^kP^^(k) + CC.) , (16) 

ak k 

£ak = 42 + ^al + ' , (17) 

~{o) _ fik^ ~(o) _ fik'' 

''^^ - 2^^' '^^-'2^^'^'' ^^^^ 
_^^y-(k-k')p"'^(k') . (19) 



^ak 



n 

k' 



The electrostatic energy eI:^1 of |$o) is zero since there is no net charge in any of the four 
layers, and the total energy ii^tot of |\E') is given by 

Etot - eSI = Ees + Y (4k + l^al) - ^ E (^kP"^(k) + CC.) . (20) 

ak ^ ^ ak 

The mean-field Hamiltonian [Eq. ( |T5l) ] is, in matrix form, 

Sm.-eS = ^,-^...pr(k)+E(4 4) f Z ■ (21) 

ak k \ ~^k ^''k / \ Ct,k / 
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For fixed pl and p^, An — rie — nh — riL + nRis fixed and we can obtain the electrostatic 
potential V^^ and the electrostatic energy £'^5 by solving a classical electrostatic problem 
of four charged layers. The electrostatic field V^^ is given by 

K'^s^^^^^ + C, (22) 

K = j^ + e,o + C, (23) 

where 

Eyo = —[PL - pR)d , (24) 

and C is a constant reference single-particle energy which we set to zero. The electrostatic 
energy per area is given up to a constant which we disregard by 

Ees 2Tie^d 



1 / 2 , 2\ fPL Pr\ 



(25) 



Now we derive the gap equation for the bilayer system with equal densities of electrons 
and holes. The system is completely defined by rii, and riR or alternatively by An and 
/^o = £vo — Eg. We use the latter pair as our parameters for the following calculations. We 
can diagonalize the Hamiltonian Hmf by introducing the new Fermi quasi-particle operators 
using a Bogoliubov transformation: 



H 



MF 




where 





(27) 



4'''^ = ^(^ck + £„k) ± ^4 + |Ak|2 , (28) 
\u^\^ = -(1+ , I , (29) 

kr = ^fi-^^^^i , (30) 



2 1 V4 + |Ak 



^^^^ ^ 2V4+|AkP ' ^''^ 

£k = (£ck — £vk) , (32) 



and the ground state is given by 

1^) = H/^klO) = n + «k4) |0) , (33) 

k k 

where |0) is the vacuum state in which both the conduction band and the valence band 
are empty. This canonical transformation is analogous to the canonical transformation of 
BCS theory if we perform an electron-hole transformation on the valence band electrons. 
In the conduction-valence band picture, these new fermion operators correspond to the 
superposition of conduction and valence band states that is induced by spontaneous phase 
coherence and diagonalizes the self-consistent mean-field Hamiltonian. In this language our 
mean-field calculation is just standard Hartree-Fock theory. 

We have described the procedure we follow in some detail because special care is required 
in how the electrostatic energy is treated in order to extract the finite energy of an electrically 
neutral system correctly. We will consistently use the conduction-valence band picture of 
electrons rather than the electron-hole picture throughout all the calculations for bilayer 
systems because it is easier to avoid confusion in how the cancelling Coulombic divergences 
are handled. The self-consistent gap equation at zero temperature is obtained from the 
definition of the inter-band exchange field, 

Ak = ^E^™(k-k')(cLc.k') 
k' 

k' 

We solved the gap equation above self-consistently for a system with an equal density 
of electrons and holes, effective masses = 0.067mo, m* = O.llmo, inter-layer distance 
d = lOOA, and dielectric constant e = 13. We choose the density difference An = He — Uh 
and the initial band overlap /io = Svo — Eg as two controllable parameters and set An = 
in this section. There is close analogy to the BCS theory of a superconductor and we obtain 
a similar energy spectrum with a gap [Figj2] . The gap function Ak has a maximum around 
the Fermi wave vector and decreases as density increases [Fig 13]. At very high densities, 
the ground state converges toward the electron-hole plasma state and the gap approaches 
zero. These results agree with the previous work by Zhu et al.^^ Even though the attractive 
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FIG. 3: (Color online) Ak in the 2-band model for various values. The calculated densities 
are n = 0.8 x lO^^cm"^, 1.0 x lO^^cm"^, 1.4 x lO^^cm"^ and 1.7 x lO^^cm"^ for = 5.7meV to 
23meV, respectively. 
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interaction is now a Coulomb interaction rather than the BCS-hke short range interaction, 
the gap function shows similar behavior, which implies that the pairing is most significant 
for electrons and holes near the Fermi surface. 

III. EXCITON CONDENSATE FERROMAGNETS 

We now restore the spin degree of freedom in the bilayer systems to study its effect 
on excitonic condensation. In population polarized systems this additional degree of free- 
dom plays an important role by allowing spontaneous spin-polarization to compensate for 
imperfect nesting of the paramagnetic state Fermi surfaces. Recently, the experimental ob- 
servation of weak ferromagnetism in lightly doped divalent hexaborides^ gave rise to new 
interest in exciton condensate ferromagnetism as a possible explanation for the experimen- 
tal finding.— 1^1^ Ferromagnetism in exciton condensates is driven by a preference for Fermi 
surface nesting of at least one spin-component. The spin-polarized state has some unique 
symmetry properties related to the approximate spin-rotational invariance of the micro- 
scopic Hamiltonian. When all operators in the theory of the previous section are assigned 
an additional spin index cr, the total Hamiltonian is 

ao-k 

+ ^ XI ^(q)cLkclv'k'Ca'<x'k'+qCaak-q • (35) 

kk'q 

aa'aa' 

As in the previous section, we define |$o) as the full- valence-band empty-conduction-band 
reference state. Electron and hole states are then specified by the density matrix 

fCU^) = (^|cLkCa.k|v&)-pSl'(k) 

= (^IcL'kCaokl^) - ($0|cL,kCa<xk|$0), (36) 

and the mean-field Hamiltonian is 

Hmf - = Eo + J2 + + ^) (4.c... - .(k)) 

acrk 
acrtr'k 

-J](Ar'cLkWk + /^.c.) , (37) 

k 
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where the intra-band exchange field has been separated into spin-independent and spin- 
dependent parts: 



^ik = E ^""(k - k')p::(k') , (38) 

o-k' 

= E ^""(k - ^'WA^')ra'a (39) 

where r is a vector whose components are Pauh spin matrices. We allow electron-hole 
pairing (conduction-valence coherence) between any spin states: 

= ^E^'"(k-k')p^;,(k'). (40) 

k' 

The band energies and the electrostatic fields are then given by the same expressions as in 
the spinless case 

V-^^-'-^ + e.o, (44) 

and the constant term is 

airk 

E^^iS-^^lk) - 1 ^ (h„k • T^.') p"A(k) 

acrk acro-'k 

+l5](Ar'p:?,(k)+c.c.) . (45) 



1 

2 



cro-'k 

The mean-field Hamiltonian in matrix form is 



Hmf = ( 4k clik citk ) H 



Ccik 
Ct,tk 
V C^ik / 



const. , (46) 
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FIG. 4: Quasi-particlc excitation energy for an 4-band cxcitonic condensate state with //q = 25meV 
and An = 0. The calculated densities are rig = n/j = 2.04 x lO^^cm"^. 



where the Hamiltonian matrix has the general form 



H 



v 



-A*r 



A *ii 

^ k 



*^k 
■ 



with 



^ - ?<0) , T/ES , r(o) 
£ak - + ^^a + "ak 



(47) 



(48) 



The total energy of ^ is given by 

+^ E (h--^ • ^-') -^vlk) - ^ E (^ry.v(k) + c.c) . (49) 

otTir'k (Tir'k 

Unlike the spinless case, the diagonalization of this 4x4 matrix is not in general triv- 
ial. In our calculations we diagonalizc it numerically and use the resulting eigenvectors 
to solve the generalized Hartree-Fock equation of the system self-consistently, allowing for 
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both spontaneous coherence and spontaneous spin-sphtting in both conduction and valence 
bands. As in the previous section, we use /xq = Svo — Eg and An = Ue—Uh as two parameters 
for the numerical calculations. When ng = rih, the self-consistent calculation converges to 
the results of the previous section except each quasiparticle state has an additional two-fold 
spin-degeneracy. (The effective magnetic field due to the intra-band exchange term h^k 
vanishes in this case.) Figure H] shows the eigenvalues of the Hartree-Fock matrix which 
are the quasi-particle excitation energies. Each quasi-particle band is doubly degenerate as 
expected. 

When He 7^ Uh the normal state Fermi surfaces of the conduction and valence bands have 
different Fermi wave vectors, frustrating inter-band coherence which is favored by bands 
with identical Fermi surfaces as we have explained. Fermi surface matching can be restored 
for one conduction-valence band pair by splitting the spin-degenerate bands so that one 
conduction-valence band pair forms a condensed state and the other pair remains in the 
normal state. This spin-polarization, however, causes the kinetic energy to increase, which 
competes with the energy gained by condensation. We therefore expect ground states with 
spontaneous coherence and spontaneous spin-polarization, provided that \ne — nh\ is not too 
small or too large. We refer to these states as exciton condensate ferromagnets. See Fig. [5] 
for a schematic illustration of this point. 

Figure [6] shows typical quasi-particle excitation energy spectrum for a system with n^. ^ 
Hh] these results are for fio = 25meV and An = 5.0 x lO^^cm^^. Spontaneous splitting 
of the spin degeneracy indeed occurs. One pair of bands show an energy gap due to the 
condensation and the other pair show no such energy gap, which indicates they have no 
interlayer coherence. Since the valence band effective mass is larger than the conduction band 
effective mass, the splitting in the valence band is much larger than that of the conduction 
band to minimize the kinetic energy cost of achieving nesting. The ground state energy and 
the energy bands are invariant under independent spin rotations in either layer. 

In bulk samples, there is a small overlap between electron and hole wave functions which is 
usually neglected (dominant term approximation). When this small overlap is not neglected, 
it gives rise to nonzero electron-hole exchange interactions, which favor spin triplet electron- 
hole pair states^ and the system is invariant only under simultaneous spin-rotations. The 
exciton condensate ferromagnet then behaves like a conventional ferromagnet with negli- 
gible magnetic anisotropy. In bilayer systems, the overlap is exponentially small and we 
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(a) 



(b) 



FIG. 5: (Color online) Schematic diagram of the Fermi surfaces and the energy bands for spon- 
taneous spin splitting in ferromagnetic exciton condensates, (a) For different number of electrons 
and holes, the conduction band Fermi surface(solid blue circle) does not coincide with the valence 
band Fermi surface (dotted red circle), (b) Conduction band electrons with spin Xc fhp to spin 
state Xc so the Fermi surface of Xc electrons shrinks and the Fermi surface of Xc electrons expands. 
Similarly, Xv valence band electrons flip to Xv state so that the Fermi surface of the Xv electrons 
increases until it matches the Xc Fermi surface. The Reconduction band and Xv valence band elec- 
trons (dot-dashed violet circle in the middle) then condense to form excitonic condensates while 
the Xc and Xv electrons remain in the normal state. The spin repopulation necessary to achieve 
Fermi surface nesting leads to ferromagnetism, i.e. to spontaneous spin polarization. 

can therefore safely neglect electron-hole exchange for typical inter-layer distances and bar- 
riers between the two layers. Since we neglect the small overlap the interaction conserves 
the band indices and the system has the spin-rotational symmetry for each band. This 
SU(2)xSU(2) symmetry leads to degeneracy of spin singlet and triplet states. The terms 
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in the Hamiltonian that break this continuous symmetry are extremely small compared to 
the bulk case. Thus we have a family of ground states which differ only in the direction 
of the spontaneous magnetization of each layer. This family of ground states has the same 
magnitude of magnetization for each layer but total magnetization magnitudes which vary 
widely. The intra-layer interaction h^^ plays an important role here, helping to stabilize 
the ferromagnetic state since it favors spin polarized states. Although the spin-polarization 
within each layer is driven by interlayer interactions, these do not create a preference for the 
relative spin orientation of the two layers. 

In numerical calculations the final converged spin-polarization direction depends on the 
initial guess we choose to begin the self-consistent loop, and different initial states lead 
to different spin-polarization directions. To verify that the pairing occurs only between one 
conduction band and one valence band leaving the other bands without coherence, we change 
the spin basis states from spin-up and spin-down to Xc and Xc for the conduction bands and 
Xv and Xv ^oi the valence bands, which are along the spin-polarizations of the conduction 
and valence bands respectively. Let the spin-orientation for the conduction (valence) band 
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FIG. 7: (Color online) The magnitudes of in different spin bases, (a) is in the spin up and 
down basis and (b) is in the new basis where the spin quantization direction is chosen to be parallel 
to the calculated total magnetization of each layer. In the new basis, the only non-vanishing A is 
l^XcXv _ This calculation is for the same model parameters used in Fig. O 

have polar angle 9c (6^) and azimuthal angle (pc {(pv)- Then the unitary operator Ul that 
connects the two spin basis sets is, 

I = f/t ( I , (50) 
where 

/ cos% sin%-e"''^'^ \ 
\ sin 1^ - cos l^e"*'^" / 
for a=c,v. The order parameters in the new basis are 

k' 
k' 

= (f^i)....^r"(f/.).„„„ , (52) 

where is Xa or Xa and aa is spin up or down. Hence we obtain. 



^k 
^k 



^XcX. ^ 








^XcX. 1 









. (53) 
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As shown in Fig. [71 only A^'^'^" is nonzero and all the other order parameters are zero in the 
new basis, which verifies that only Xc and Xv spin bands are coherent. 

The Hamiltonian can now be separated into a normal part and a condensed part: 



Hmf = Hn + Hc , (54) 



where 



\. J • '''' 



and h^^ is the Xa-parallel component of the intra-layer exchange field in the new spin basis, 
which is the only non-vanishing component of the exchange field. The condensate part can 
be diagonalized using the Bogoliubov transformation as in the 2-band model studied in the 
previous section, which leads to the eigenstates 

|lk) = Icxck) , (57) 

|2k) = \vxv^ , (58) 

|3k) = Mklcxck) - vl\vxv^) , (59) 

|4k) = Wklcxck) + Mk|t^Xt,k) , (60) 

with eigenvalues 

£ik = £ck - h^k , (61) 

£2k = e.k - hJk , (62) 

£3k = , (63) 

£4k = . (64) 



where 



^'■'^ = \ (^ck + e.k + h,^^ + hJi:) ± ^4+|AkP , (65) 
e^. = \ (^ck - ^.k + h^k - Kl) ^ (66) 
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and 




(67) 



(69) 



(68) 



The ground state is given by 



n 4 n 411410) 



(70) 




where |0) is the vacuum state with no electrons and kpc (krv) is the Fermi wavevector for 
the normal conduction (valence) band. 

When the density difference is very small, the pair of bands with different Fermi surfaces 
can also form a condensate and the spin repopulation readjusts to minimize the total energy. 
In the other limit, if the density difference is very large, the kinetic energy cost accompanying 
the spin repopulation will be larger than the condensation energy gain. Then there will be 
no condensation at all and the system will remain in the normal paramagnetic state. The 
states discussed above in which only one condensate occurs, are found in a broad intermediate 
regime. The states which manage this compromise between different energy contributions 
are found automatically by solving the mean-field equations. 

IV. INFLUENCE OF RASHBA SPIN-ORBIT INTERACTIONS 

In 2D layers, the Rashba SO interaction appears due to the structural inversion asymme- 
try of the confining potential. The strength of the SO coupling can be tuned by applying an 
external electric field perpendicular to the layers.— '^'^'^ The Rashba Hamiltonian can be 
derived using Lowdin perturbation theory^i^ up to the third order.— The effective Rashba 
SO interactions for the conduction band and for the heavy hole valence band are respectively 




ak X z ■ T 




(71) 
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and 



■ r , (72) 



P 



where r's are the Pauh matrices, t± = (r^. ± iTy)/2, k± = kx ± iky, and tan0k = ky/kx- 
These SO interactions can be considered as momentum- dependent effective magnetic fields 

hfk = a{ky^-kxy) , (73) 
hi = (sin 30kX - cos 30ky) , (74) 

which change direction in spin-space as directions changes in momentum space [see Figs. [H^a) 
and[9]^a)]. For the conduction band, the energy dispersion gets an additional linear term 
e'^^ = ± ak and the corresponding spin states are 

|ck±) = -^(|ckt)T^e^'^1cki)) , (75) 

which are shown in Fig. MJo). For the heavy hole valence band, we get an additional term 
proportional to k^, e''^ = ± (3k'^ and the corresponding spin states are 

\vk±) = i= {\vk t) T te''^-\vk i)) , (76) 

which are shown in Fig. Mjo). 

Unlike the case of the electron gas without Rashba SO interactions, the ground state of 
the free electron gas does not solve the Hartree-Fock equation of the system with Rashba 
SO interaction. The spin states for each k are determined by the total effective magnetic 
field h^^ which is the sum of the intra-band exchange field h^k and the Rashba SO field h^. 
For the higher energy band with energy = Eaii + jh^kl? ^Pi^ parallel to the total 
effective magnetic field while for the lower energy band with energy el~^ = Saii — |h^f |, the 
spin state is antiparallel to the total effective magnetic field. The self-consistency conditions 
for the intra-band exchange field then leads to 

~ 1 

h-'c = E ^ (k - k') {Pl\ - ir.) , (77) 
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(a) (b) 

FIG. 8: (Color online) (a) Rashba SO effective magnetic field and (b) the spin states for 
the conduction band. The spin direction is (p^ — vr/2 for |ck+) (inner circle and blue arrows) and 
0k + 7'"/2 for |ck—) (outer circle and red arrows). 




(a) (b) 

FIG. 9: (Color online) (a) Rashba SO effective magnetic field and (b) the spin states for the 
heavy hole valence band. The spin direction is 3(/>k — 7r/2 for |wk+) (outer circle and blue arrows) 
and 3(/)k + 7r/2 for |t>k—) (inner circle and red arrows). 
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where the spinor of the state |ak±) is parallel(+) or antiparallel(-) to the total effective field 
h^f , not just the Rashba field. The total spin of band a is given by 

1 

= 2 5Z ™ (Pr+(k') - P-".(k')) . (78) 

In the normal state the effect of exchange interactions is simply to enhance the magni- 
tude of the Rashba interaction induced spin-splitting of the bands. For exciton-condensate 
ferromagnets, however, the Rashba and exchange fields are not in general parallel. 

The Rashba SO Hamiltonian breaks spin-rotational invariance around an arbitrary axis, 
but maintains invariance under simultaneous spin and orbital rotations around the z- 
direction, and under inversion z —>■ —z combined with inversion in the xy plane (k — k). 
Thus we expect that the energy of the exciton condensate ferromagnet is independent of the 
azimuthal angle 0a of the magnetization, and invariant under reversal of the polar projec- 
tion, i.e. Qa — Qa- When the in-plane components of the Rasbha and exchange fields are 
not parallel, quasiparticle energies will depend on momentum space orientation. There will 
therefore tend to be some orientations in which the conduction and valence band Fermi en- 
ergies are close and some orientations where they are more widely separated. Condensation 
then occurs mainly in the region in which the two Fermi surfaces are close together. 

In the following we show that, depending on densities and on the strength of the Rashba 
SO interaction, the total spin of the ground state can be along the growth direction (uniaxial 
ferromagnet) or the total spin can have nonzero in-plane component. When the total spin is 
along the z direction, we calculate the anisotropic energy by applying an external magnetic 
field to force the total spin direction off the z axis and subtracting the magnetization energy. 
When the total spin has a nonzero in-plane component, the total spin of the ground state can 
have any azimuthal angle. The rotational symmetry around the z direction is spontaneously 
broken in this case and the quasiparticle energy spectrum is not rotationally symmetric in 
k space. We show representative examples for each case below. 

Experimental gate control should enable experimental control over the relative strength 
of the Rashba interactions in valence and conduction bands.— In order to reveal the basic 
physics more simply we consider systems in which the Rashba SO interaction acts only in 
one layer. Figure [10] shows our results when the total spin is along the z direction. The left 
panels in this figure are for a system with SO interaction only in the conduction band. The 
parameters are a = 0.05 eVA, (3 = 0, hq = 30meV, and An = Ue — nh = 4.0 x lO^^cm^^. The 
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densities obtained self-consistently are = 2.59 x lO^^cm"^ and Uh = 2.19 x lO^^cm"^. The 
direction of the total spin of the conduction band layer is 6c = 0.047r and (pc = 1.727r, which is 
almost along the z direction. In this uniaxial case, the quasiparticle spectrum is rotationally 
symmetric in momentum space [Fig. [10] (b) and (c)]. This circumstance is achieved in the 
self-consistent solution by converging to a state in which the intra-band exchange field for 
the conduction band hck has in-plane components which are parallel to the Rashba field h^. 
The total effective field of the conduction band has the form 

= h^{k) (sin0kx - cos</.ky) + K{k)z . (79) 

This form implies that the magnitude of the total effective magnetic field is independent 
of the momentum space orientation (/)k, explaining the isotropic Fermi surfaces. In these 
solutions, the in-plane component of the total spin [Eq. (ITHI) ] for the conduction band layer 
vanishes and we get total magnetization Sc only along the z-direction. On the other hand, 
the absence of SO interaction in the valence band leads to arbitrary direction for the effective 
exchange field and therefore the total spin Sy for the valence band layer is in an arbitrary 
direction. Similar considerations apply for a = and (3 [right panels in Fig. [10] , where 
we use a = 0, P = 500 eVA^, /iq = 30meV, and An = rie - Uh = -4.0 x lO^^cm^^. The 
calculated densities are ng = 2.24 x lO^^cm"^ and Uh = 2.64 x lO^^cm"^. The direction of 
the total spin of the valence band layer is 6^ = 0.997r and 0^, = 0.857r. The self-consistent 
solutions converge to states with concentric constant-energy surfaces and the total effective 
field is of the form 

hll = hx(fc) (sin30kX - cos30ky) + Kik)z . (80) 

Again, the in-plane component of the total spin [Eq. (1781) ] for the valence band layer vanishes 
and we get total magnetization along the 2;-direction. The direction of the total spin of the 
conduction band layer is arbitrary due to the absence of SO interaction in the conduction 
band. 

For these uniaxial cases, we can calculate the total energy of the system with magnetiza- 
tion directions off the ^-direction by introducing weak external magnetic field. Specifically, 
we apply an external magnetic field strong enough so that magnetization is off the z axis by 
da, then subtract the magnetization energy M ■ Hext from the total energy to calculate the 
anisotropic energy of this uniaxial ferromagnet. Figure [TT] shows the numerically evaluated 
total energy for (a)Q: 7^ 0, /5 = and (b)a = 0, /9 7^ 0, as a function of cos^ 9a. It shows 
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that the total energy of the uniaxial system is roughly proportional to cos^ 9a- The energy 
difference between the state with magnetization aligned to the z axis and the state with 
magnetization perpendicular to the z axis is a few /ieV. Because the magnetic anisotropy 
energy is very small, Fig. [TT] looks somewhat noisy. We can also derive this linear behavior 
by treating the Rashba SO interaction as a perturbation (see Appendix [X] for details). The 
zeroth order ground state is given by Eq. (!70|) . We calculate the perturbed energy for each 
quasi-particle state using the states Eq. fl57j) ~ Eq. fl60|) as the unperturbed quasi-particle 
states. The total perturbed energy is evaluated by summing the corrections for each quasi- 
particle state up to the Fermi energy of the unperturbed ground state, assuming that the 
Fermi energy does not change much by the Rashba SO interaction. The first order correction 
vanishes. At second order we obtain 

Etot = + Aa^ cos^ Oc + 5/5^ cos^ 9^ , (81) 

where 6E is the energy correction that does not depend on the magnetization angles and A, 
B are constants. 

Even though a perturbative calculation indicates that the ferromagnet should have either 
an easy axis or an easy plane, our non-perturbative self-consistent calculations sometimes 
find ground state with neither an easy axis nor easy plane. Instead, in some cases the 
ferromagnet can have a non-trivial optimal polar angle and an arbitrary azimuthal angle. 
Evidently higher order contributions can change the ground state qualitatively. For a given 
spontaneously chosen azimuthal angle the quasi-particle band structure is anisotropic. This 
broken XY symmetry leads to intricate and quite interesting quasi-particle properties. The 
left panels of Fig. [12] show a case with a 7^ and /5 = where the total spin of the 
conduction band has a nonzero in-plane component. The parameters used for this system 
are a = 0.03 eVA, (3 = 0, fio = 30meV, and An = He - nn = 4.0 x lO^^cm'^. The 
calculated densities are ng = 2.59 x lO^^cm^^ and Uh = 2.19 x lO^^cm^^. The direction of 
the total spin of the conduction band layer is 6^ = 0.457r and (pc = 1.727r. As can be seen 
in Fig. [TW b) and (c), the conduction band quasi-particle excitation energy dispersions are 
not rotationally symmetric. The constant-energy surfaces for the conduction band shift so 
that the two bands are closer in one direction and farther apart in the opposite direction in 
k space. This corresponds to an intra-band field that points in the same direction as the 
Rashba field, but the magnitude of the intra-band field is a function of not only of k but 
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also of 0k- To get nonzero in-plane components of the total spin along the 0c direction, 
the total effective magnetic field is stronger when the azimuthal angle of the spin states of 
the majority species(|ck— ) for the conduction band) is 0c and weaker when the azimuthal 
angle of the spin states of the minority species(|ck+) for the conduction band) is 0c. The 
constant-energy surfaces are farther in the direction where the effective field is stronger and 
closer where the effective field is weaker since the energy difference between the two spin 
bands is the magnitude of the effective magnetic field. Therefore, we obtain the condition 
for the closer (farther) constant-energy surfaces 



The blue dashed arrows in Fig. [T2|b) and (c) show the 0k that satisfy the closer condition 
above, which agrees with the numerically calculated constant-energy surfaces. 

The right hand panels of Fig. [12] show a case with a = and /3 7^ for which the 
valence band magnetization has a nonzero in-plane component . The parameters used for this 
calculation were (3 = 400 eVA^, a = 0, /iq = 20meV, and An = — Uh = —4.0 x 10^°cm~^. 
The calculated densities are rif. = 1.59 x lO^^cm"^ and = 1.99 x lO^^cm^^. The direction 
of the total spin of the valence band layer is 6^ = O.GStt and 0^ = O.SStt. For the valence 
band, the majority species has spin state |fk— ) and minority species has |f k+) because the 
valence band has opposite energy dispersion curve compared to the conduction band. So 
the majority band has spin states that point to 30k + 7r/2 and the minority species has spin 
states that point to 30k — n/2. Thus there are 3 closer directions and 3 farther directions 
that satisfy 



Fig. [T2r e) and (f ) show the directions of the closer condition by blue dashed arrows, which 
again agrees with the numerical results very well. These systems spontaneously break the ro- 
tational symmetry around the ^-direction to make the Fermi surfaces of one conduction band 
and one valence band come close together so that they can form the excitonic condensate 
while the other two bands remain normal. 



0c = 0k - f ; closer 
0c = 0k + f ; farther 



(82) 



0^ + 2mT = 30k — f ; closer 
(py + 2mi = 30k + I ; farther 



(83) 
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V. SUMMARY AND DISCUSSION 



We have studied electron-hole pair condensation in spatially separated bilayer systems. 
Our work is motivated by recent progress^'^^'^'^ in the preparation of electron-hole bilayers 
with carriers which are generated electrically rather than optically, and in equilibrium rather 
than in a steady state. These systems already show behavior, in particular enhanced drag 
voltages at low-temperatures, which appears to be evidence for a non-Fermi liquid ground 
state driven by attractive interactions between conduction band electrons and valence band 
holes. The extremely large drag voltages expected^ in the two-dimensional superfluid state 
have not yet been seen, suggesting that the Kosterlitz-Thouless temperatures of current 
samples are still below available temperatures. More robust experimental consequences can 
be anticipated if systems can be fabricated in which electron-hole interactions have been 
strengthened by reducing the quantum well widths and hence the average distance between 
layers. 

The most unique and attractive aspects of equilibrium electrically generated electron-hole 
systems are i) the opportunity to directly probe how transport properties are altered^ by 
excitonic superfluidity and ii) the opportunity to continuously adjust the relative density of 
electrons and holes. The study of pairing fermion systems with unbalanced populations^ 
has been a major topic in cold atom physics. In this paper we have explored some of the 
peculiarities expected to be associated with unbalanced populations in the case of electron- 
hole bilayers. The main difference between cold-atom systems and electron-hole bilayers 
single-particle states have an additional attached spin-label in each layer. This additional 
degree of freedom is expected to qualitatively change the population polarization physics. 
Population polarization and the subsequent Fermi surface mismatch induces spontaneous 
spin-polarization which improves nesting between some Fermi surface pieces and therefore 
increases the pairing condensation energy. When spin-orbit interactions are neglected, the 
spin-polarized state is invariant under independent spin-rotations conduction valence band 
layers. In this paper we have examined how Rashba SO interactions, which will inevitably 
be present in any equilibrium quantum well bilayer system, alter these ferromagnetic con- 
densate states. The Rashba SO interaction arises naturally because of the electric field 
used to transfer electrons between layers. The interaction introduces a structural inversion 
asymmetry which breaks the spin-rotational symmetry and leads to a preferred polar angle 
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of the magnetization of each layer. When an in-plane component of the magnetization is 
present it implies anisotropic energy spectra and interaction driven spatial anisotropy in 
the two-dimensional plane which should be readily detected experimentally. Because the 
anisotropy energy is very small compared to other energy scales, and also because of small 
crystalline anisotropy effects which are not included in our calculation and appear when the 
spin-orbit interaction is expanded to higher order in two-dimensional momentum, we have 
not been able to establish patterns in the relationship between magnetization direction and 
the experimentally controllable system parameters like the density, the density difference 
between electrons and holes, and the strength of the Rashba SO interaction. 

Work at the University of Texas was supported by the National Science Foundation under 
grant DMR-0606489. 



APPENDIX A: PERTURBATION THEORY FOR EXCITON CONDENSATE 
FERROMANGET WITH RASHBA SO INTERACTION 

The mean-field Hamiltonian of the ferromagnetic excitonic condensate of electron-hole 
bilayer system without SO interaction is given by 



H 



MF 



Hn + Hi 



c 



(Al) 



wherei^AT is the Hamiltonian for the normal components [Eq. (155|) ]. and He is for the con- 
densate [Eq. fl56|) ]. For each k, we have 4 quasiparticle eigenstates given by Eqs. fl57|) ~ 
( I60l) . with eigenvalues Eqs. (I6T1) ~ ( I64l) . These quasiparticle eigenstates are related to the 
conduction and valence band states by 



/ Cctk \ 



I 



smf 

cos f e^<^ 



= U 






C2k 
C3k 
V C4k / 







sin ^■ 

cos %e* 



cos y 

sin l^e*"^ 



f k cos ^ 
-vl^ sin \e 



f k COS Y \ [ Clk \ 

f k sin 1^6*'^'= C2k 

Uk cos ^ C3k 

Mk sin l^e'"^" / \ C4k / 



(A2) 



(A3) 
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The ground state of this system is given by Eq. (ITOll . We will consider this state as our 
unperturbed state and treat the Rashba SO interaction as a perturbation. 

In the basis of {|c t k), |c | k), \v ] k), \v [ k)} Rashba SO Hamiltonian is given by 



H 



R 



/ iake-'^ 
-iake'"^^ 

\ 

and, in new basis {|lk), |2k), |3k), |4k)}, by 











i(3k^e-^''^^ 
-i(3k^e^'^^ 



Hr' = U+HrU . 



(A4) 



(A5) 



Applying perturbation theory for each k, first order terms for the quasiparticle energies 



are 



-(1) 

'Ik 

-(1) 
'2k 

-(1) 
'3k 

-(1) 
'4k 



Hr j = ak sin 9^ sin(</)c - </)k) , 
Hr' ) = [3k^ sin 6^ sin(0i, - 30k) , 

/ 22 

= —ak \u\S^ sin 6c sin(0c — 0k) — l3k^ \v\S^ sin 0^ sin( 
= —ak Itikp sin Oc sin(0c — 0k) — l3k^ |-Uk|^ sin 0^ sin( 



H 



R 
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H 



R 
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and the change in the total energy is 

4 



5E, 



(1) 

tot 



'ik / ^ik 



EE® i^^-^ 

1=1 k 



(A6) 
(A7) 
30k) , (A8) 
30k) , (A9) 



(AlO) 



where the integration over 0k vanishes for all e-i"*. So there is no contribution from the first 



order terms. Second order energy corrections for each k are calculated by 



.(2) 



fH 



JO) _ JO) ■ 



(All) 
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After some algebra, we get 
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2 ; 2 





+ 








£31 




£41 
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£42 



[cos^(0c - 0k) + cos^ 9c sin^(0c - 0k)] , 
^ [cos^(0^ - 30k) + cos^ 9y sin^(0^ - 30k)] 



Z 7 z I I z 

(2) Ot rv l^k r 2// /\ 2/1 - 2// /\1 
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where Sij — ef^ — ef^ Using 



/ (i0k sin^(0c - 0k) = / #k cos^(0c - 0k) = tt , 

JQ Jo 
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we can find second order energy correction 
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(A19) 



Spin-rotational symmetry is now broken and depending on the sign of the constants A and 
The total spin can have either easy plane or easy axis. 
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FIG. 10: (Color online) Exciton condensate Ferromagnet with uniaxial magnetization with Rashba 
SO interaction in one of the layers. Left panels [(a)~ (c)] are for a = 0.05 eVA and f3 = 0. The 
direction of the total spin of the conduction band layer is along the z direction, (a) shows the 
quasiparticle energy spectrum and the dashed black line is the chemical potential, (b) and (c) are 
constant-energy surfaces in momentum space, corresponding to the energies shown in (a), (b) is 
for 0.025 eV (red dot-dashed line in (a)) and (c) is for 0.012 eV (blue dotted line in (a)). Right 
panels [(d) ~ (f)] are for /5 = 500 eVA^ and a = 0. The direction of the total spin of the valence 
band layer is along the z direction, (d) shows the quasiparticle energy spectrum (e) and (f) are 
constant-energy surfaces in momentum space, (e) is for 0.023 eV (red dot-dashed line in (d)) and 
(f) is for 0.01 eV (blue dotted line in (d)). 
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FIG. 11: Magnetic anisotropy of uniaxial systems. /Uo=30meV and An=5.0 x lO^'^cm"^. (a)a = 
0.05 eVA , P = 0, (h)a = 0, (3 = 700 eVA^. With these parameters, the ground state magnetization 
is along the z direction. We apply small external magnetic field to change the direction of the 
magnetization and then subtract the magnetization energy to extract the energy of the system. 
The energy per particle shows a linear dependence on the cos^ 9a where 6a is the angle between 
the z axis and the total spin for band a = c,v. 
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FIG. 12: (Color online) Exciton condensate ferromagnet with nonzero in-plane component of the 
magnetization with Rashba SO interaction in one of the layers. Left panels [(a)~(c)] are for 
a = 0.03 eVA and /? = 0, where the direction of the total spin of the conduction band layer 
is 9c = 0.457r and (j)c = 1.727r.. (a) shows the quasiparticle energy spectrum and the dashed 
black line is the chemical potential, (b) and (c) are constant-energy surfaces in momentum space, 
corresponding to the energies shown in (a), (b) is for 0.025 eV (red dot-dashed line in (a)) and (c) 
is for 0.014 eV (blue dotted line in (a)). The black solid arrow depicts the direction of the total 
spin in the conduction layer and the blue dashed arrow depicts the direction of i?!>k = + 7r/2. 
Right panels [(d)~(f)] are for [i = 400 eVA^ and a = 0, where the direction of the total spin of 
the valence band layer is 6^ = O.GSvr and (j)^ = O.SSvr. (e) and (f) are constant-energy surfaces in 
momentum space, (e) is for 0.015 eV (red dot-dashed line in (d)) and (f) is for 0.005 eV (blue 
dotted line in (d)). The black solid arrow depicts the direction of the total spin in the valence band 
layer and the blue dashed arrows depict the direction 0k = 0t,/3 + 7r/6 + 2n7r/3 for integer n. 
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